Universal three-body physics for fermionic dipoles 



O 

O 

CO 

3. 

Ph. 



O 

C/5 ■ 

Ph. 



> 

m 
m 

o 



Yujun Wang, J. P. D'Incao, and Chris H. Greene 
Department of Physics and JILA, University of Colorado, Boulder, Colorado, 80309-0440, USA 

A study of the universal physics for three oriented fermionic dipoles in the hyperspherical adiabatic 
representation predicts a single long-lived three-dipole state, which exists in only one three-body 
symmetry, should form near a two-dipole resonance. Our analysis reveals the spatial configuration of 
the universal state, and the scaling of its binding energy and lifetime with the strength of the dipolar 
interaction. In addition, three-body recombination of fermionic dipoles is found to be important 
even at ultracold energies. An additional finding is that an effective long-range repulsion arises 
between a dipole and a dipolar dimer that is tunable via dipolar interactions. 
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The recent advances in producing ground-state ultra- 
cold dipolar molecules [1] have sparked great interest in 
the novel phases that could be experimentally accessi- 
ble [2]. With the ability to control the dipolar inter- 
action by an external electric field [3], studies of dipo- 
lar molecules could stimulate the discovery of intriguing 
phenomena across a broad range of disciplines, from con- 
densed matter physics to ultracold chemistry. In par- 
ticular, the effects of the long-range anisotropic dipolar 
interaction on ultracold physics has already been studied 
for three bosonic dipoles [4, 5]. A recent development 
has shown [5] that the Efimov effect [6] is more favorable 
for three bosonic dipoles although the dipolar interaction 
has ingredients that could conceivably have destroyed it. 

Stimulated by their ability to mimic a variety of ex- 
otic physical systems from superfluids to neutron stars, 
quantum degenerate fcrmi gases with nondipolar inter- 
actions have been studied extensively in recent years [7] . 
In contrast to bosonic gases whose lifetime is mainly 
determined by three-body recombination [8], fermionic 
gases exhibit extraordinary stability against three-body 
inelastic collisions except when the two-body interaction 
is tuned near a p- wave resonance [9, 10]. Moreover, three- 
body physics for identical fermions is expected to be 
nonuniversal, i.e., short-range-dependent [9, 11]. Con- 
sequently, three-body physics for identical fermions have 
to some extent been regarded as less exciting and/or dif- 
ficult to observe. Nevertheless, the influence of dipolar 
interactions on the universal behavior of bosons triggers 
new interest in the study of universality for fermionic 
dipoles, and it is becoming increasingly important to un- 
derstand such physics in order to guide ongoing experi- 
ments on ground-state fermionic dipolar molecules [3]. 

In this Letter, the behavior of three identical fermionic 
dipoles is demonstrated to be universal, thus differing 
fundamentally from the nondipolar cases. In particular, 
near a dipole-dipole resonance (i.e., where the two-dipole 
binding energy E^d — > 0), a single long-lived, universal 
three-dipole state is formed. Its size and lifetime both 
grow with stronger dipolar interactions. The stability of 
the fermionic dipolar gases is assessed through numeri- 
cal calculations of the rate coefficient K3 for three-dipole 



recombination, D + D + D D2 + D, which is the dom- 
inant expected loss mechanism for nonreactive ground- 
state dipolar molecules [12]. In fact, we have found that 
K3 grows rapidly with the strength of the dipolar interac- 
tions, and it is strongly enhanced near a dipole-dipole res- 
onance. This study also shows evidence for a long-range, 
dipolar interaction-dependent repulsion that suppresses 
D-2 + D relaxation collisions when D2 is deeply-bound. 

Consider the three-dipole Schrodinger equation in hy- 
perspherical coordinates (in atomic units): 
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H ad Up = Eip. 



(1) 



where ip = i? 5 / 2 ^ is the scaled three-body wavefunction, 
/x=m/v3 is the three-body reduced mass for identical 
particles with mass m. The adiabatic H a d is 
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where A 2 is the usual grand angular momentum opera- 
tor [13], and Q are the hyperangles representing positions 
of the three particles at fixed hyperradius R. The poten- 
tial energy V = v(fi2) + v(r23) + v{r23) is a pairwise sum 
of two-dipole interactions: 



v(f) = Vbsech 2 (r/7- ) + 
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f(r). (3) 



The sech term is isotropic and short-range, and it en- 
capsulates the complicated interactions between two po- 
lar molecules. The second term with a short-range cutoff 
f(r) is the interaction between two dipoles aligned in the 
z direction. The dipole length is defined in terms of the 
electric dipole moment d m as di = md 2 n /2. The univer- 
sality of our calculations is tested by varying both di and 
V . 

Equation (1) can be routinely solved, either using the 
adiabatic representation [13] or else through slow variable 
discretization (SVD) [14]. Our study shows it is advanta- 
geous to use the adiabatic scheme to treat the asymptotic 
behavior of three dipoles in scattering calculations, while 
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SVD efficiently handles sharp avoided crossings at small 
distances. In both cases, the biggest challenge is to solve 
the adiabatic eigenvalue equation: 



U v {R)^ M {R-n) 



(4) 



where U V (R) is the adiabatic potential for channel u, 
M is the space-fixed frame projection of the total angu- 
lar momentum J, and &^ ,M (R;tt) is the corresponding 
eigenfunction. For interacting oriented dipoles, the ma- 
jor difficulty in solving Eq. (4) is that the total orbital 
angular momentum J is not conserved. The method of 
Ref. [5] is implemented to solve Eq. (4). Briefly, $JJ ,M is 
expanded in terms of the Wigner D functions [see Eq. (4) 
in Ref. [5]], truncated at J max = 13. Tests of this trun- 
cation confirm that it yields adequate convergence. 

The nature of three-dipole states near a dipole-dipole 
resonances emerges from study of the adiabatic poten- 
tials for AI n = + . Note that the even and the odd J 
values decouple in this symmetry. This symmetry with 
only odd J is of particular interest because it includes the 
J u = 1 + , the least repulsive partial wave for three nonin- 
teracting identical fcrmions [16]. One would also expect 
the binding of three dipoles is most likely to occur at 
mid = dipole-dipole resonance, where m-id is the pro- 
jection of two-dipole angular momentum along the field 
direction. However, rri2d = is not allowed for A/ n = + 
symmetry, so the three-dipole system is considered in- 
stead near the vn-id = 1 resonances. Figure 1 shows the 
typical behavior of the adiabatic potentials. Asymptoti- 
cally, the three-body continuum potentials behave in the 
same manner as the nondipolar interaction case, but for 
dipole plus dipolar dimer channels, each threshold has 
a family of potential curves whose centrifugal barriers 
carry different effective angular momenta due to the J- 
coupling. 

Irregardless of the complicated topology of the adia- 
batic potentials shown in Fig. 1 , near a dipole-dipole res- 
onance it is possible to trace a diabatic potential well in 
the channel closest to the three-body continuum. It cuts 
through multiple sharp avoided crossings with deeply- 
bound channels, as is visible in the inset of Fig. 1. Re- 
markably, all these potentials for different dg and for dif- 
ferent short-range potentials fall on top of each other 
after a proper scaling. A key implication of the scal- 
ing of the universal potentials is that the hyperradius of 
the repulsive barrier increases with dg, which suppresses 
decay of the three-dipole state and increases its lifetime 
as the dipolar interaction is increased. To quantify the 
number of three-dipole states and their energies in the 
universal potential wells, we have solved Eq. (1) by using 
SVD. The energy E 3 d and the width T of the three-body 
states are revealed by adding complex absorbing poten- 
tials in the lower decaying channels. Interestingly, ex- 
actly one quasi-bound state is supported by each of the 
universal potentials. Table I lists E$d and V for a few 
values of dt. The energy E^d shown in Table I suggests 




FIG. 1: (color online) A typical set of adiabatic hyper- 
spherical potentials U V (R) for three fermionic dipoles with 
dg/rO — 58.2 and E 2 d — > 0. The inset shows rescaled, diaba- 
tized potentials exhibiting universal behavior for a few values 
of de at a dipole-dipole resonance. The horizontal dashed line 
in the diabatic potential wells indicates the position of the 
universal three-dipole states. 
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TABLE I: The upper rows: the energies E^d and the widths 
r of the universal three-dipole states for different values of 
dt. The bottom row: the angle between the three-dipole 
triangle and field direction, the long and short bond length bi 
and 6 S , and the smaller bond angle A from the most probable 
configuration of the three-dipole states. 



a universal trend, but the width T exhibits less universal 
behavior, presumably due to the nonuniversal couplings 
to deeply-bound channels. Nevertheless, the sup- 
pression in the width implies an increased lifetime of the 
universal three-dipole states as dg increases. Our numeri- 
cal study also shows that the universal three-dipole states 
slowly become unbound or deeply-bound as the system 
shifts away from a dipole-dipole resonance. In a broad 
region near a dipole-dipole resonance, the properties of 
the three-dipole states are essentially unchanged. 

Importantly, the interactions of the aligned dipoles 
break the overall rotational symmetry, reflecting that 
the three-dipole states have a preferential orientation in 
space. The geometry of these universal states is analyzed 
by calculating the spatial distribution of the dipoles using 
the three-body wavefunctions. Specifically, the probabil- 
ity density for finding a dipole at r, is 



p(?) = \(*\J26(r-r i )\*), 



(5) 
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FIG. 2: (color online) The upper part: the isosurfaces of 
the one-dipole probability density p(f) for the universal three- 
dipole states, centered in the space-fixed center of mass frame. 
By normalizing the peak density to be unity, the density of 
the following surfaces: the red (innermost), three green (mid- 
dle), and the light blue (outermost) are 0.8, 0.3 and 0.015, 
respectively. This crankshaft-resembling structure is a figure 
of rotation about the field axis. The lower part: the most 
probable orientation and geometry of the triangle formed by 
three dipoles. The green axes show the direction of the exter- 
nal electric field. 



where 5(f) is the Dirac delta function, f*, is the position 
vector for the i'th dipole. Figure 2 shows the isosurfaces 
of p(r). The three-dipole states have a strongly prefer- 
ential spatial configuration, as indicated by the triangle 
in the lower part of the figure. The isosurfaces can be 
viewed as the revolution of the triangle along the field 
axis. It is also interesting to note that the configuration 
of the three-dipole states are universal, only their overall 
size scales linearly with dg. 

The configuration of a three-dipole state shown in 
Fig. 2 can be understood in terms of the anisotropic na- 
ture of the dipolar interactions and the symmetry of the 
state. The three dipoles tend to align in a line where the 
dipolar interactions in Eq. (3) are the most attractive, 
but this is forbidden by the odd partial wave symmetry. 
The state therefore has the configuration of an obtuse, 
isosceles triangle making a small angle with the field di- 
rection. The lower part of Table I lists the angle between 
the plane of the triangle and the field direction, the short 
and long bond lengths, and the smaller bond angle de- 
termined from the most probable and also the dominant 
configuration of the three-dipole states. 

The possibility of universal three-dipole states for other 
symmetries have been explored, in particular the three- 
dipole system for M n = 0~ symmetry with odd J's, 
which includes the less repulsive partial wave contribu- 
tion J n = 1~ [16]. However, no universal three-dipole 
states have been found near both m.2d = and m,2d = 1 
dipole-dipole resonances. The situation is similar for 
AI n = 1~ and 1 + symmetries, where the adiabatic po- 



tentials are found to be more repulsive. The potentials 
for higher M near higher mid resonances are expected 
to be more repulsive, suggesting the conclusion that the 
M n = + three-dipole states presented above are the 
only class of universal states for identical fermions. 

The universal three-dipole states affect scattering pro- 
cesses. As dg increases, a resonance in three-dipole re- 
combination is expected whenever a three-dipole state 
passes through the three-body breakup threshold before 
a new two-dipole state becomes bound. To determine 
the stability of nonreactive fermionic dipolar gases, how- 
ever, it is more important to know the general scaling 
of three-dipole recombination with dg and the behavior 
near a dipole-dipole resonance. The threshold behavior 
of the recombination rate K 3 for three fermionic dipoles 
can be obtained from a Wigner law analysis [16] of the 
AI U = + , J-odd low-energy-dominant contribution, giv- 
ing 

K 3 = (6) 



where k = \J1[iE is the three-body wavenumber. 

For short-range, nondipolar interactions, the coeffi- 
cient C3 in recombination rate is expected to scale with 
the scattering volume [9, 10] and the effective range [10]. 
For dipolar interactions, however, the scattering physics 
is fundamentally different even in the two-body level [17]. 
In particular, we have numerically identified a distinct 
low-energy expansion of the p-wave phase shift 5 for two 
fermionic dipoles. The real part of 6 can be expanded as 

Re[<5(fc 2 )] ~ ~ak 2 - bk\ - Vk\ (k 2 -> 0), (7) 

where fc 2 is the relative two-dipole wavenumber. This 
expansion has a linear term in k 2 implying a nonvanishing 
scattering length a [17] as expected for r -3 interaction 
and also a quadratic dependence on k 2 , both of which 
are absent for nondipolar interactions. But similar to 
the nondipolar interaction case, the scattering volume V 
goes through a pole across a dipole-dipole resonance. 

Our scaling laws for three-dipole recombination are 
based on a calculation of the three-dipole scattering 
matrix using the eigenchannel i?-matrix propagation 
method [18]. First, SVD solution of Eq. (1) obtains the 
i?-matrix out to a distance that includes all the sharp- 
avoided crossings. The i?-matrix is then propagated far- 
ther in the adiabatic representation [13] in the asymptotic 
region until K3 is fully converged. K3 has been calcu- 
lated for V > near a dipole-dipole resonance, where 
the dipolar dimer is weakly bound. At ultracold ener- 
gies, our calculation shows that only the lowest three- 
body continuum channel is important for recombination, 
and K3 is dominated by recombination into the weakly- 
bound dipole plus dipolar dimer channel having the low- 
est effective angular momentum barrier asymptotically. 
Figure 3 shows the scaling coefficient C3 near mid = 1 
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FIG. 3: (color online) The dependence of the recombination 
scaling coefficient Ca on the two-body scattering volume V 
and the dipole length de near a dipole-dipole resonance. The 
lines connecting the numerical points just guide the eye. 



dipole-dipole resonances for a few di with different short- 
range interactions. The dipole-dipole resonances with 
higher m,2d are not of current interest because of their 
extremely narrow widths [17]. Dimensional analysis has 
suggested, and our numerical results in Fig. 3 confirm, 
that a universal scaling of C3 exists 



C 3 = \V 17/2 d, 



-35/2 



(8) 



where dimensionlcss A « 2 x 10 9 . This scaling is stronger 
with V than for nondipolar interactions [9, 10]. Although 
Eq. (8) also seems to imply a strong suppression of three- 
dipole recombination with di, instead a general scaling 
of C3 oc d^ is expected away from a dipole-dipole reso- 
nance, due to the d\ scaling of the background value of 
V [19]. This rapid growth of C3 with di may be unfortu- 
nate for experiments that require stable fcrmionic dipolar 
gas in three dimensions, since three-dipole recombination 
can be dramatic even in the ultracold threshold regime 
(kdi <C 1). On the other hand, scenarios exist where 
three-body recombination has a positive role, as a means 
to produce weakly-bound dipolar dimers that are stable 
against collision with a free dipole, in close analogy to 
the technique used in ultracold spin-mixed fermions with 
nondipolar interactions [20]. 

Finally, consider the behavior of deeply-bound dipole 
plus dipolar dimer channels and their impact on the sta- 
bility of a mixed gas containing both dipolcs and dipolar 
dimers. Similar to the bosonic dipole case [5], these chan- 
nels have a centrifugal barrier characterized by a non- 
integer effective angular momentum quantum number, 
and it grows with increasing de. Although our numeri- 
cal study shows that this (^-dependence is nonuniversal, 
the effective angular momentum barrier can nevertheless 
"protect" a mixture of dipoles and dipolar dimers from 
collisional decay, and it appears to provide an alternative 
way to tune interactions in a mixed-species gas. 

In summary, a class of universal physics has been iden- 
tified for three fcrmionic dipolcs. For scattering volumes 



near a dipole-dipole resonance, one universal thrce-dipolc 
state exists at an energy near the three-dipole breakup 
threshold. This state has a universal geometry and orien- 
tation in the external field, with a variable size dependent 
on di. Moreover, its long lifetime can be beneficial for 
experimental observation and manipulation. Our study 
also shows a universal trend of the three-dipole recom- 
bination rate, which increases rapidly near each dipole- 
dipole resonance, and grows rapidly with stronger dipolar 
interactions. And the effective repulsion between a dipole 
and a dipolar dimer can be tuned to control collisions in 
mixed quantum gases. 
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